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In this review, I present a brief summary of the impact of nucleon pairing 
at supra-nuclear densities on the cooling of neutron stars. I also describe 
how the recent observation of the cooling of the neutron star in the 
supernova remnant Cassiopeia A may provide us with the first direct 
evidence for the occurrence of such pairing. It also implies a size of the 
neutron 3 P — i-2 energy gap of the order of 0.1 MeV. 

1. Introduction 

With radii of the order of 10 km and masses between 1 and 2 M©, neutron 
stars a giant nuclei the size of a city. Neutron star matter at densities 
around /9 n uc, the nuclear saturation density, is certainly made of neutrons 
and protons, with electron and muons to guarantee charge neutrality. At 
higher densities it is likely that hyperons and/or deconfined quarks are 
present, as well as charged meson condensates. The presence of such "ex- 
otica" can have a dramatic effect on the equation of state (EOS), softening 
it and reducing the maximum possible neutron star mass. The recent mea- 
surement of a two-solar-mass neutron star is, however, imposing very severe 
constraints on the EOS [1] . The chemical composition of dense matter also 
has a strong effect on the cooling of a neutron star, by determining the 
allowed neutrino emission channels, a fact that provides another handle on 
the EOS. However, the occurrence of pairing, which strongly affects exci- 
tations but has only a very weak effect on the EOS, dramatically alters the 
neutrino emissivity of the star and hence its cooling. Comparison of the- 
oretical models with observations almost certainly requires the occurrence 
of pairing, and the recent observation of the rapid cooling of the neutron 
star in Cassiopeia A [2] may be the first direct evidence for it. 
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In this chapter, I will present a brief description of the physical processes 
relevant for the modeling of neutron star cooling followed by a series of 
generic results, and an overview of the "Minimal Cooling" paradigm and 
its application to Cas A. This presentation is necessarily short and the 
reader can find more details in the reviews [3j-[5]. An introduction to the 
physics of neutron stars can be found in [5] . 



2. The Basic Physics of Neutron Star Cooling 

Most neutron star cooling calculations involve solving numerically the heat 
transport and energy balance equations, in their fully general-relativistic 
form (see, e.g., Appendix B in [7]). However, simple global energy-balance 
(in its Newtonian form to keep it simple) as 

dEfh „ dT , . 

-ir = Cv^ = -L v -L 1 + H (1) 

illustrates the main features, where it is assumed that the star's interior 
is isothermal with temperature T, a state reached within a few decades 
after its birth in a core-collapse supernova. Here, is the star's total 
thermal energy, Cy its specific heat, and L 7 and L v its photon and neutrino 
luminosities, respectively. The term H , for "heating" , represents possible 
dissipative processes which will not be considered here. 



2.1. Specific heat 

The core provides most of the specific heat and since it consists of degen- 
erate fermions one naturally has Cy oc T. When matter is made only of 
nucleons and leptons, about 70% of Cy is provided by the neutrons, 20% 
by the protons and 10% by the leptons. For simple numerical estimates 
one can use Cy ~ 10 39 Tg erg K _1 , where Tg = T/10 9 K. Notice, however, 
that in the presence of pairing Cy can be strongly suppressed. 

2.2. Photon emission and the envelope 

The surface photon luminosity is traditionally written as 

L 7 = 4ttR 2 a SB T e 4 (2) 

which defines the effective temperature T e (R is the star's radius and ctsb 
the Stefan-Boltzmann constant). T e is much lower than the star's interior 
temperature: even in the case the interior is isothermal there always exits 
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Tabic 1. Dominant neutrino emission processes. 
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t The coefficients i?'s are control functions to incorporate the effects of pairing, see § 12,41 

a strong temperature gradient in the uppermost layers, commonly called 
the envelope. As a rule of thumb one can use 

T e w 10 6 (T/10 8 K) 1/2 K (3) 

which implies that L 7 « 10 35 Tg erg s _1 . The details of this T e — T relation- 
ship depend on the chemical composition of the envelope and the presence 
or absence of a strong magnetic field. 

2.3. Neutrino emission processes 

A list of the most important neutrino emission processes is presented in 
Table [TJ with rough values of their emissivities. They are separated into 
"slow" and "fast" processes, the former involving 5 and the latter only 3 
degenerate fermions. Notice the different temperature dependences: T 6 for 
the three degenerate fermion processes compared to T 8 for the five fermion 
ones, a direct consequence of the stronger phase space limitation resulting 
in a significantly reduced emissivity. The Cooper pair process is described 
in § 12.41 A detailed description of neutrino processes can be found in 8 
and an alternative approach in [9]. 
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2.4. Effects of nucleon pairing 

The occurrence of pairing, either of neutron or of protons 0, introduces a 
series of important effects: 

A) Alteration, and possible strong suppression when T <C T c , of the specific 

heat of the paired component. 

B) Reduction, and possible strong suppression when T <C T c , of the emis- 

sivity of neutrino processes the paired component is involved in. 

C) Triggering of the "Cooper pair breaking and formation" (PBF) neutrino 

process which is very efficient in the case of spin-triplet pairing. 

These effects are direct consequences of the alteration of the quasi-particle 
spectrum by the development of an energy gap D(k), 

normal phase : e(fc) = vp(k — kp) (4) 

paired phase : e(fe) = ±y / [vp(k — kp)] 2 + D 2 (kp) (5) 

which severely limits the available phase space when T <Si T c . In cooling 
calculations these effects are introduced through "control functions" : 

cy — > R c cy and e v — > Rx ^ ■ (6) 

There is a large family of such control functions for the various types of 
pairing and the numerous neutrino processes "X" (they are the factors 
R in Table [IJ. For nodeless gaps the R's are Boltzmann-like factors ~ 
exp[— 2A(T)/fcsT] and result in a strong suppression when T <C T c , while 
for gaps with nodes the suppression is much milder. Regarding the specific 
heat, there is a sudden increase, by a factor ~ 2.4 at T = T c , followed by a 
reduction at lower T. 

The effect C), neutrino emission from the formation and breaking of 
Cooper pairs |T0l Qj], can be seen as an inter-band transition (as, e.g., 
n — > n + vv) where a neutron/proton quasiparticle from the upper (+) 
branch of the spectrum of Eq. <j5j> falls into a hole in the lower (— ) branch 
Such a reaction is kinematically forbidden by the excitation spectrum of the 
normal phase, Eq. but becomes possible in the presence of an energy- 
gap, Eq. (JSJ) - The resulting emissivity can be significantly larger than the 
one of the modified Urea process in the case of spin-triplet pairing. 



a Hyperons, and deconfined quarks, if present, are also expected to pair. I will, here, only 
consider nucleons but most of what follows naturally translates to these others cases. 
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2.5. Theoretical predictions of nucleon gaps 

Pairing is usually assumed to occur in a single spin-angular momentum 
channel A = (s,j) and the gap function is a 2 x 2 matrix in spin space A\. 
At low k Fn , or k Fp , it is theoretically predicted that the preferred channel 
is A = (0,0) in S-w&ve, i.e., the spin-singlet Sq. At larger Fermi momenta 
the So interaction becomes repulsive and the preferred channel is A = (1, 2) 
in P and F waves (the mixing being due to the tensor interaction), i.e., the 
spin-triplet 3 P— F2. In the So channel, which has also been called the "A" 
phase, the gap is spherically symmetric and Amo) depends on one single 
scalar A(fc) so that the energy gap D(k F ), at the Fermi surface, is simply 

A phase (^o) : D{k F ) = A(k F ) . (7) 

In the 3 P— F2 channel, A\ has contributions from all possible orbital angular 
momentum I and rrij components, i.e., A\ — m , A™£ (fc)G^ (k F ) where 
the (k F ) are 2x2 spin matrices describing the angular dependence of A 
which is thus not spherically symmetric. Microscopic calculations restricted 
to the 3 P2 channel indicate that the largest component of A\ corresponds 
to the rrij = sub-channel or, possibly, the rrij = ±2 one, sometimes called 
the "B" and "C" phases, respectively. For these two special cases the energy 
gap D(kp) is given by [12] 

B phase ( 3 P 2 ,m, = 0) : D 2 {k F ) = \ [A» A (fc F )f 1 + ° (8) 

C phase ( 3 P 2 , mj ; = ±2) : D 2 (k F ) = [A 2 2X (k F )] 2 (9) 

where 9 is the angle between kp and the arbitrary quantization axis. The 
relationship between the phase transition critical temperature T c and the 
energy gap D(k F ) is approximately given by the usual result 

k B T c w 0.57 A(fc F ;T = 0) (10) 

for all three phases [12j H3j . where A(k F ;T = 0) is obtained by angle 
averaging of D 2 (k F ) over the Fermi surface 

[A(k F ;T = 0f^ JJ^D 2 (k F ;T = 0). (11) 

Figure [T] shows a sample of theoretical predictions of T c for the neutron 
3 P — F2 gap. The three dotted lines present some of the first published 
models: "HGRR" from TF], "T" from [T5] and "AO" from [12 . The four 
continuous lines show models from |19j and illustrate the uncertainty on 
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Fig. 1. Left panel: some theoretical predictions of T c for the neutron i P — F2 gap. See 
text for description. Right panel: some phenomenological models of T c for the neutron 
3 P — F2 gap used in neutron star cooling simulations. Models "a" , "b" , and "c" are 
from [7] and |14| . model "a2" from |15l . On the top margin are marked the values of 
k Fn at the center of a 1.0, 1.4, 1.8, and 2.0 M star built with the APR EOS [16]. 



the gap size due to the problem that no N — N potential reproduces the 
measured 3 P2 phase-shift above E\ a t> ~ 300 MeV (translating into fc^ n ~1.8 
fm _1 ). These gaps were calculated with the Nijmegen II ("NijH"), Ni- 
jmegen I ("NijI"), CD-Bonn ("CDB"), and Argonne Vis ("AVI 8") poten- 
tials (displayed values are taken from the middle panel of Figure 4 of [1.9| ) . 
These calculations are at the "BCS approximation" level, i.e., do not in- 
clude medium polarization effects. In the case of the gap, medium 
polarization is known to result in screening and to reduce the size of the 
gap. In the case of a 3 P2 gap, polarization with central forces is expected 
to result in anti-screening and to increase the size of the gap. However, 
Schwenk & Friman [20j showed that spin-dependent non-central forces do 
the opposite and strongly screen the coupling in the 3 P2 channel, resulting 
in a T c lower than 10 7 K: this "SF" value is indicated in the Figure by an 
arrow ! Finally, 3-body forces are know to be essential for both nuclear 
structure and neutron star matter. They are repulsive in the bulk but at 
the Fermi surface in the 3 P — channel they turn out to be strongly at- 
tractive. The two dashed lines in Figure [T] present results from [3T| where 
the "2bf" model only considers 2-body forces (from the Argonne Vis) while 
the "3bf" includes a meson exchange model 3-body force: the result is a 
growing 3 P— gap which shows no tendency to saturate at high density. 
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Fig. 2. Some theoretical predictions of T c for the proton 1 So gap. See [7] for references. 
On the top margin are marked the values of kp p at the center of a 1.0, 1.4, 1.8, and 2.0 
M© star built with the APR EOS [16]. 

A set of theoretical predictions of T c for the proton x So gap is shown in 
FigureO Variations in the size of this gap are not as large as for the neutron 
triplet gap but uncertainty on the kp range is significant, translating in 
an uncertainty of a factor ~ 4 on the density range covered by proton 
superconductivity. 

Given these large uncertainties on the size of the neutron 3 P — F% gap 
(about three orders of magnitude) and the fact that neutrino suppression 
depends on it through an exponential Boltzmann-like factor, this gap is 
often considered as a free parameter in neutron star cooling models. The 
extreme sensitivity of the cooling history on the size of this gap can be 
taken to one's advantage, by inverting the problem, as it may allows us 
to measure it by fitting theoretical models to observational data [22] • The 
right panel of Figure Q] presents the phenomenological neutron 3 P—F 2 gaps 
used in the following cooling calculations. 

3. Some Illustrative Examples of Neutron Star Cooling 

To illustrate several of the possible behaviors of a cooling neutron star, I 
show in Figure [3] a set of models, all based on the same star of mass 1.4M Q 
built with the APR EOS (TBJ- The "slow cooling" models only include, in 
the core, the slow neutrino processes of Table[T]as well as the BPF process. 
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Fig. 3. Left panel: cooling curves for various illustrative cooling scenarios. Right panel: 
temperature profile evolution for the fast cooling model with q = 26. The numbers on 
the curves give the age of the star, in years. See text for description. 



For the "fast cooling" models I added, at p > 3p nuc , a fast process with 
emissivity e^ st = 10 q ■ T 9 6 erg cm- 3 s _1 , with q = 25, 26, and 27, that 
simulates the effect of a kaon condensate, a pion condensate, or a direct 
Urea, respectively. These models, all based on the same EOS, are not self- 
consistent but have the advantage that the only differences between them 
is the presence/absence of the e^ st process and the presence/absence of 
pairing. The models with pairing include the neutron 5*0 gap from |20] in 
the inner crust, a ^So proton gap in the outer core from [23], model "T" of 
Figure [2j and the phenomenological neutron 3 P — F2 gap "b" of Figure [T] 

The various distinctive phases of the evolution are marked as "A" , "B" , 
"C", and "D" above the cooling curves. Phases A and B are determined 
by the evolution of the curst while C and D reflect the evolution of core. 
Phase A: the surface temperature T e is determined by the evolution of the 
outer crust only. At such early stages the temperature profile in the outer 
crust is independent of what is happening deeper in the star and all models 
have the same T e . 

Phase B: the age of the star becomes comparable to the thermal relaxation 
time-scale of the crust and heat flow controls the evolution of T e . This 
thermal relaxation phase is depicted in the right panel of the Figure which 
shows the evolution of the T-profile for the fast cooling model with q = 26 in 



March 8, 2013 13:36 World Scientific Review Volume - 9in x 6in Pairing'Cooling 



Pairing and the Cooling of Neutron Stars 9 

the absence of pairing (marked as "Normal" in the left panel). One sees that 
very early on the "pit" (where the fast neutrino emission is occurring) is very 
cold and during the first 30 years heat is flowing from the outer core into 
the pit whose temperature consequently remains stationary. Afterward, 
the core is essentially isothermal and heat from the curst is now rapidly 
flowing into the cold core, a process which takes little more than 100 years 
and is reflected by a sudden drop of T e when the cooling wave reaches the 
surface. After this, during phases C and D, the stellar interior is isothermal 
and it is only within the shallow envelope, not shown in the Figure, that a 
temperature gradient is still present. 

Notice that the models with pairing have a shorter crust relaxation time: 
this is simply due to the strong reduction of the neutron specific heat in 
the inner crust by the x So gap. 

Phase C: during this phase, the "neutrino cooling phase", the star evo- 
lution is driven by neutrino emission from the core since L v 3> L 7 . The 
difference between "slow" and "fast" neutrino emission is now clearly seen. 

Very noticeable is the effect of the pairing suppression of the neutrino 
emissivity in the fast cooling models: once T dropped below T c , which 
happened only a few seconds/minutes after the beginning of the simulation, 
neutrino cooling was quenched. Notice that the difference between the 
various e~p ast processes is much smaller in the presence of pairing than in 
its absence: that it took half a minute (for q = 27) or half an hour (for 
q = 25) for T in the "pit" to drop below T c does not make much of a 
difference when looking at the star thousands of years later. The evolution 
is more dependent on T c than on the actual value of q. 

Considering the slow cooling models, comparison of the normal case 
with the paired one, but with the PBF processed turned off, shows the 
same effect of pairing, but less dramatic since the star is not as cold as in 
the fast cooling models. However, when the PBF process is not artificially 
turned off the paired model is very similar to the normal one: the burst 
of neutrino emission occurring when T ~ T c from the constant formation 
and breaking of Cooper pairs induces some extra, temporary, cooling. The 
impact of the PBF process, however, depends on the size of the neutron 
P—F2 gap and is considered in more detail in the next section. 
Phase D: at late times, L v has dropped significantly due to its strong T 
dependence and photon emission, L 1 , now drives the evolution. This is 
seen in the cooling curves by their very large slopes. During this "photon 
cooling era" the models with pairing cool faster due to the reduction of the 
specific heat. 
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4. The Minimal Cooling Paradigm and Cas A 

Table [1] shows that there are many possibilities for having some fast neu- 
trino emission and, from the results of the previous section, this implies 
that any observed T e could be fit if there is a gap of the appropriate size 
that controls the neutrino emission. However, a simple question arise: is 
there any observational evidence for the occurrence of fast neutrino cool- 
ing? If yes, how strong is it? To address this question, the "Minimal 
Cooling" paradigm was developed [3 [14] which is a natural extension of 
the previous "Standard Cooling" (see, e.g., [21] )• The essence of the mini- 
mal cooling paradigm is the a priori exclusion of all fast neutrino emission 
processes. Core neutrino cooling is, hence, limited to the modified Urea 
with the similar nucleon bremsstrahlung processes, and the PBF process. 
Minimal cooling is, however, not naive cooling: it takes into account all 
other uncertainties on the microphysics and the astrophysical conditions. 
The major factors turn out to be the size of the neutron triplet-pairing gap 
and the chemical composition of the envelope. 

It was shown in Figure[3]that, with a relatively large value of the neutron 
3 P—F2 gap, the PBF process can compensate the suppression of the modified 
Urea process. A further study of this effect is presented in Figure [U In 
the left panel, a cooling curve with a vanishing 3 P — F2 gap and three 
curves with the gaps "a", "b", and "c" of Figure [1] are presented. The all- 
important parameter is the maximum value of T c reached in the core, T™ ax , 
since the pairing phase transition will start when the core temperature 
T reaches T™ ax . It is seen that colder stars during the neutrino cooling 



T 




Age [yrs] Age [yrs] 

Fig. 4. Gauging the effect of neutrino emission by Cooper pair formation (the PBF 
process) within minimal cooling. Large values, > 10 9 K, of for the neutron 3 P— i*2 

gap are considered in the left panel while smaller values, < 10 9 K, are considered in the 
right panel. See text for description. 
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era are obtained with a value T™ 3 * = 10 9 K than for larger values. The 
reason for this result is that with such a gap as model "a" there is always 
a significant region of the core that is going through the phase transition 
at ages between 10 to 10 5 years while in the cases of the larger gaps "b" or 
"c" the pairing phase transition occurred within most of the core at much 
earlier times and most of the neutrino emission, including the one from the 
PBF process, is hence suppressed at later ages, resulting in warmer stars 
during the neutrino cooling era. The right panel of Figure [4] describes the 
effect of smaller neutron 3 P — F2 gaps: they are the model gap "a" scaled 
by a factor x = 0.6, 0.4, or 0.2, resulting in T™ ax = x x 10 9 K. Ones sees 
that the cooling trajectories for x < 1 separate from the upper trajectory, 
with no 3 P - F 2 gap, at a point which is precisely the moment when the 
core temperature T reaches T™ ax and the pairing phase transition starts El- 
This late onset of the pairing phase transition results in a transitory period 
of accelerated cooling. 

4.1. Minimal cooling vs data 

The present set of observational data on isolated cooling neutron stars is 
displayed in the Figures [5] and [B] The first subset of stars, displayed as 
boxes and numbered 1 to 13, are objects from which a thermal spectrum, 
in the soft X-ray band Ex ^0.1 — 5 keV, is clearly detected. For each 
case the plotted T e is obtained from the thermal luminosity using Eq. ((2J. 
The stars are: 1) CXO J232327.8+584842 (Cas A), 2) PSR J1119-6127 
(G292. 2-0.5), 3) PSR J0821-4300 (Puppis A), 4) PSR 1E1207.4-5209 (PKS 
1209-52), 5) PSR B0833-45 (Vela), 6) PSR B1706-44 (G343. 1-2.3), 7) PSR 
B0538+2817 (S147), 8) PSR B2334+61(G114.3+0.3), 9) PSR B0656+14, 
10) PSR B1055-52, 11) PSR B0633+1748 aka "Geminga", 12) RX J1856.5- 
3754, and 13) RX J0720. 4-3125 (in parenthesis is given the name of the 
associated supernova remnant). The second subset consists of four pulsars 
whose soft X-ray spectrum is dominated by magnetospheric emission, ap- 
pearing as a power-law spectrum, so that only an upper limit on T e can 
be inferred. They are: A) PSR B0531+21 (Crab), B) PSR J1124-5916 
(G292.0+1.8), C) PSR J0205+6449 (3C58), and D) PSR J0007.0+7303 
(CTAl). Finally, six upper limits, shown as dotted arrows, are obtained 
from the total absence of detection of any object in six supernova remnants: 
some of these remnants may contain an isolated black hole, but it is un- 



b A model with x = 0.8 would be essentially undistinguishable from the case x = 1 as 
the phase transition would start before the end of the crust relaxation phase. 
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Neutron 3 P-F 2 gap: "a" (T C T = 10 9 K) Neutron 3 P-F 2 gap: "b" (T™ x =3xl0 9 K) 




Age [yrs] Age [yrs] 

Fig. 5. Comparison of two minimal cooling scenarios with observational data. The 
neutron 3 P — Fi gaps employed are shown in Figure [T] See text for description. 



likely to be the case for the six of them. I refer the reader to [7] and [14] 
for details and discussion of these data points. 

In Figure [5] these data are compared with two series of models with a 
moderate (left panel) and a large (right panel) neutron 3 P — F 2 gap. The 
grey shaded areas, encompassing several cooling curves, show the range of 
predictions resulting from the uncertainty on the chemical composition of 
the envelope: the curve marked "Fe" assumes the presence of iron peak 
nuclei while the other curves show the effect of having increasing amounts 
of light elements in the envelope (the interior physics being the same for all 
curves). Stars with ages less than 10 5 years are the most critical to infer 
properties of dense matter since they are in the neutrino cooling era. In 
the case Tj^ ax = 10 9 K, the theoretical predictions are compatible with all 
observed values and upper limits, with the only exception of the PSR in 
CTAl (D) if one does not consider the six non-detection from supernova 
remnants. On the other hand, for larger T™ ax theoretical predictions are 
incompatible with the majority of data points during the neutrino cooling 
era. 

Given the uncertainty on the theoretical predictions of the size of the 
neutron 3 P— i<2 gap, these results are not very encouraging: a moderate gap 
implies that most observed young neutron stars show no evidence of fast 
neutrino emission while a large gap results in the opposite conclusion that 
most observed young neutron stars show evidence of fast neutrino emission. 

The older stars, 10 to 13, clearly appear warmer that the theoretical 
predictions: they likely imply the presence of some "heating" mechanism, 
the "H" term in Eq. (JTJ (see, e.g., [4]). 



March 8, 2013 13:36 World Scientific Review Volume - 9in x 6in Pairing'Cooling 



Pairing and the Cooling of Neutron Stars 13 

4.2. The cooling of the neutron star in Cassiopeia A 

To be able to distinguish between the various possible neutron 3 P—F 2 gaps 
some new observational datum(a) is needed. This has likely been recently 
provided by the observation of the cooling, in real time, of the neutron star 
in the Cassiopeia A supernova remnant [2] ■ Between the year 2000 till 2009 
this neutron star appears to have cooled by 4% (with a 20% decrease in 
flux). The observed thermal luminosity of Cas A is L 1 ~ 10 34 erg s _1 : from 
the observed cooling rate, and a simple estimate of its specific heat, Cas A 
must be loosing about 10 38 erg s^ 1 , in the case the cooling represents 
the evolution of an isothermal star. It is hard to imagine that such an 
enormous cooling rate is a persistent evolution, and it is more likely that 
"something critical" happened recently: for a cooling neutron star this 
"something" should be the temperature and a critical temperature means 
a phase transition. The right panel of Figure |4] precisely exhibits such 
sudden increase in the cooling rate when the star's core temperature reaches 
T c ™ ax . With the known age of Cas A of 330 years B, the right panel of 
Figure S] indicates a T™ ax of the order of 0.5 x 10 9 K. This interpretation 
was proposed in [15] and, simultaneously and independently, in [25]. The 
left panel of Figure [6] is an up-date of Figure [5] using the neutron 3 P — F 2 
gap model "a2" of Figure Q] as proposed in [15] . 

The deduced value of T c ™ ax is mostly fixed just by the known age of the 
star. Reproducing the observed cooling slope, s = ~d\ogT e /d\ogt ~ 1.2, 
however, requires the star was quite hot before the onset of the phase tran- 
sition, which means its previous neutrino luminosity was low. Suppression 
of L v is naturally obtained if protons were already superconducting, with 
a higher T c . Obtaining s ~ 1.2 requires proton superconductivity in a 
significant part of the core, and this observed slope would also put strong 
constraints on proton superconductivity: results of Figure [5] assume the 
proton gap "CCDK" of Figure [5J Within the set of proton ^So gaps 
presented in Figure[5]this model is the one reaching to the highest densities. 
The first report of [2] was subsequently consolidated by another observa- 
tion reported in [25) . However, there are delicate calibration issues [26) and 
the actual slope may be smaller, a possibility that would ease the strong 
requirement on the proton So gap and/or allow for a more massive star. 
The right panel of Figure [6] shows the change in the evolution when increas- 
ing the neutron star mass: a larger mass simply results in having a larger 
portion of the core where protons are not superconducting and, hence, a 



c The supernova explosion was likely observed by J. Flamsteed on 16 August 1680. 



March 8, 2013 13:36 



World Scientific Review Volume - 9in x 6in 



Pairing'Cooling 



ii 



Dany Page 



Neutron 3 P-F 2 gap: "a2" (T c max = 5.5x10 S K) 




1 10 100 1000 10 4 10 5 10 6 10 7 
Age [yrs] 





M/M 9 T C [10'K] 
-1.9 0.51 : 
—1.6 0.52 ~ 
pl.3 0.57 ■ 





200 500 
Age [yrs] 



Fig. 6. Left panel: a minimal cooling scenario reproducing the observed rapid cooling 
of Cas A (data point 1). The inset, from |15| . shows the fit to the five temperatures 
reported in [2]. Right panel: dependence of the slope s of the cooling curve on the star 
mass: s = 0.5, 0.9, and 1.4 for M = 1.3, 1.6, and 1.9A/ Q , resp., at 330 yrs (from [15]). 

higher initial L v and a smaller slope during the rapid cooling phase. Sim- 
ilar results are obtained by changing the density range covered by proton 
superconductivity while keeping the star mass constant. 

5. Conclusions 



If the cooling of Cas A is confirmed by future observations, it is certainly 
one of the most amazing pieces of observational data on neutron stars and 
the interpretation presented here would mean that we are seeing, in real 
time, about a solar mass of neutrons going through the triplet pairing phase 
transition. It would also imply a neutron P — F2 energy gap of the order 
of 0.1 MeV. 

In a more general context, almost any cooling scenario compatible with 
observational data needs some type of pairing, either for protons, neutrons 
or exotica. A possible exception being the "medium-modified Urea" sce- 
nario [9]. Given the large number of possible scenarios it is difficult, when 
going beyond the minimal cooling paradigm, to decide which type of ex- 
otica is acting (see, e.g., [27J ). However, some cold young neutron stars, 
as the one in CTA1, tell us that more than minimal cooling is definitely 
needed. 
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